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We investigate a possibility to generate non-classical states in light-matter coupled noisy quan¬ 
tum systems, namely the anisotropic Rabi and Dicke models. In these hybrid quantum systems a 
competing influence of coherent internal dynamics and environment induced dissipation drives the 
system into non-equilibrium steady states (NESSs). Explicitly, for the anisotropic Rabi model the 
steady state is given by an incoherent mixture of two states of opposite parities, but as each parity 
state displays light-matter entanglement we also find that the full state is entangled. Eurthermore, 
as a natural extension of the anisotropic Rabi model to an infinite spin subsystem, we next explored 
the NESS of the anisotropic Dicke model. The NESS of this linearized Dicke model is also an in¬ 
separable state of light and matter. With an aim to enrich the dynamics beyond the sustainable 
entanglement found for the NESS of these hybrid quantum systems, we also propose to combine 
an all-optical feedback strategy for quantum state protection and for establishing quantum control 
in these systems. Our present work further elucidates the relevance of such hybrid open quantum 
systems for potential applications in quantum architectures. 
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I. INTRODUCTION 

Light-matter coupled quantum systems are now seen 
as novel composite systems for the exploration of vari¬ 
ous tasks in quantum computation and information pro¬ 
cessing, as well as for testing foundations of quantum 
physics m- These hybrid systems promise to com¬ 
bine disparate quantum degrees of freedom in construct¬ 
ing scalable quantum architectures ma- In particular, 
the idea is to take ‘the best of two worlds’: matter pro¬ 
vides good candidates for storage of quantum informa¬ 
tion, while light (photons) are superior when it comes 
to transmitting or sending quantum information, with 
interactions between these provide processing. For ex¬ 
ample, photons can be sent between different parties in 
order to create entanglement over macroscopic distances. 
However, physical quantum systems, in this case espe¬ 
cially the light part, invariably also couple to their ex¬ 
ternal environments mi. It is commonly argued that 
non-classical states, including entangled states, are ex¬ 
tremely sensitive to noise and dissipation. For exam¬ 
ple, environment induced decoherence tends to reduce 
quantum-coherent superpositions to incoherent mixtures 
[7]. Quantum state engineering strategies have therefore 
taken a center stage in salvaging quantum coherence in 
hybrid quantum systems miHi. It is clear that for using 
photons as information carriers between different mat¬ 
ter subsystems, sustainable entanglement between light 
and matter is a necessity, i.e. such quantum correlations 
should survive any realistic decoherence/noise affecting 
the photons. 
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In a somewhat parallel approach, noisy coupled quan¬ 
tum systems have been proposed as exciting avenues for 
controlled generation of multi-partite entangled states 
[9]. The interplay between coherent and incoherent dy¬ 
namics in interacting multi-partite open quantum sys¬ 
tems can result in the generation of steady states which 
exhibit exotic quantum character mm To take a par¬ 
ticular example relevant for this work, for open optical 
systems the dynamics can often be well described by a 
Lindblad type master equation 

p = -i[Hs,p]+'yC^p, ( 1 ) 

where p is the system’s state, Hs the system Hamilto¬ 
nian, and jC^p is the Lindblad super operator render¬ 
ing the effects of the reservoir on the system {A is the 
so called quantum jump operator, and we note that Hs 
may be modified from the closed system Hamiltonian due 
to the interaction with the environment). The steady 
state is solved by letting i[Hs^Pss] — = 0- It fol¬ 

lows that if pss is a dark state, i.e. C^Pss = 0, then 
Pss must also commute with Hs, e.g. being an eigen¬ 
state of the Hamiltonian. However, this is not the only 
possibility for having a steady state; both contributions 
might be non-zero but also exactly cancel each other. 
This would mean that the unitary evolution ‘balances’ 
the non-unitary evolution. For infinite systems, this in¬ 
terplay may result in non-equilibrium phase transitions 
between phases supported either by the first or the sec¬ 
ond term in Eq.([T) (TOl E]. That is to say that pss 
changes qualitatively and shows non-analytic properties 
at some critical 7c. For finite systems, on the other 
hand, such a transition is smooth/analytic, but never¬ 
theless for the intermediate stages between the two ex¬ 
tremes non-classical states may persist. In this work 
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we explore this scenario and examine the possibility of 
generating non-classical states in two open quantum sys¬ 
tems that have served as work horses in quantum optics 
for more than half a century, namely the Rabi m (or 
non-Rotating-Wave-Approximation (non-RWA) Jaynes- 
Cummings model [13]) and the Dicke model [14]. In a 
way, in our study these two models are at the two ex¬ 
tremes of the quantum spectrum] the Rabi model which 
couples a single spin-1/2 particle to a boson or photon 
mode lives in the deep quantum regime, while what we 
term the Dicke model considers instead a particle with 
an infinite large spin S (almost classical) coupled to the 
boson mode. As already mentioned, the main source 
of decoherence is very often photon absorption and we 
thereby consider a Lindblad term with A = a being the 
photon annihilation operator. This photon loss mecha¬ 
nism aims at emptying the boson mode (the only dark 
state is the vacuum). However, for a large light-matter 
coupling in the Hamiltonian even the ground state (in 
the absence of environmental coupling) becomes a ‘highly 
excited’ (from the perspective of photon number) state, 
and so the two dynamical contributions in the evolution 
of the hybrid system must be balanced in order to sup¬ 
port a NESS. We will demonstrate that as a result of this 
competition between unitary and non-unitary evolution, 
for both models a NESS exists that possess non-classical 
features in terms of light-matter entanglement. Eurther- 
more, especially for the Rabi model this NESS reflects 
the well known parity symmetry of the original Hamil¬ 
tonian; the NESS consists of an incoherent mixture of 
different parity states. 

In the second half of the paper we explore possibilities 
to systematically enrich and enhance the non-classical 
properties of the NESS. More precisely, although envi¬ 
ronment induced dissipation potentially can be tailored 
to achieve desired quantum states [15] , establishing quan¬ 
tum control over various dissipation channels is also an 
important benchmark toward realizing a useful quantum 
architecture m- Thus, it is desirable to construct quan¬ 
tum control schemes for hybrid quantum processors if, for 
instance, one is interested in initial quantum state preser¬ 
vation m- Along these lines we will use an all-optical 
coherent feedback strategy [T8H2Q] to achieve quantum 
state preservation. It is important to appreciate that 
such a scheme is different from just tailoring the dissipa¬ 
tion channels as it involves active feedback to the system. 
As outlined in Refs. mms], under the assumption that 
the time-delay introduced by the feedback loop is negli¬ 
gible it is indeed possible to establish a complete control 
over the dissipative dynamics. Going beyond such ideal¬ 
ized sudden feedback control, we will conclude by briefly 
commenting on the prospects of our scheme by combin¬ 
ing it with a time-delayed feedback control scheme as 
previously studied both in the classical [21] and quan¬ 
tum domains [22] [23] . 

The paper is organized as follows. In Secjlljwe intro¬ 
duce our first physical model describing the interaction 
between a quantized bosonic mode and a single two-level 


system. We explore the general structure of the NESS p^s 
(being of an interesting form) and extract the sustainable 
entanglement between the two subsystems. The follow¬ 
ing Sec |HI| addresses the question whether non-classical 
features survive in the ‘classical limit’ in which the two- 
level system is turned into an infinite one i.e. identifying 
the linearized Dicke model. And indeed, entanglement is 
sustained also in this case. Having discussed these two 
model examples, in Sec. [^ we combine a coherent all- 
optical feedback scheme for quantum state protection in 
open quantum systems. While entanglement is sustain¬ 
able in the original open systems, the feedback scheme 
allows for increasing the desired properties. Einally, we 
conclude the paper with a short discussion in Sec. |V| 

II. MODEL A - ANISOTROPIC RABI MODEL 

The quantum Rabi model describes the interaction be¬ 
tween a quantized bosonic field and the simplest quan¬ 
tum mechanical model, a single two-level system or qubit 
[12] . Despite its simplicity, the quantum Rabi model had 
eluded an exact solution for many decades and it was only 
recently that an analytical solution was obtained [24] . 
with continuing debate as to the completeness of this so¬ 
lution and the existence of other possible solutions [25] . 
Nonetheless, over the years the quantum Rabi model has 
turned out to be an ubiquitous physical model and has 
found applications in understanding a wide variety of 
physical systems, including trapped ions, superconduct¬ 
ing qubits, optical and microwave cavity QED, circuit 
QED, among others [26]. In some cases the two-level 
system may be an exact description of the matter part 
of the hybrid system (when it is indeed a spin-1/2 par¬ 
ticle), but in many cases the two-level description is an 
approximation of a more complicated matter sub-system. 

In its most simplified form the quantum Rabi model 
takes the form (with h=l and neglecting the vacuum en¬ 
ergy of the boson mode) 

^RM = cjd’f'd + + A(d‘*‘ + d)((j+ + a~). (2) 

Here, d and are the annihilation and creation operators 
for the bosonic field of frequency cj, = {ax ± 
with ax^y^z fhe Pauli matrices for the two-level system, Q 
is the energy level splitting between the two levels, and 
the coupling strength between the bosonic mode and the 
two-level system is represented as A (> 0). Through¬ 
out we will use dimensionless parameters such that we 
scale all frequencies by uj and time by Neverthe¬ 

less, we keep uj in all expressions in order to keep track 
of all terms and instead always use cj = 1 in any calcu¬ 
lations. The counter-rotating terms in the Hamiltonian 
([^, a^a~^ and (j“d, do not conserve the excitation num¬ 
ber. However, if the coupling A ^ cj, 11, these counter¬ 
rotating terms can be dropped from the Hamiltonian 
Q under the so-called Rotating Wave Approximation 
^WA) [131128] . Simplifying the Hamiltonian (§ under 


3 


the RWA results in the widely studied Jaynes-Cummings 
model which can be straightforwardly solved nsi i3o]. 
As evident through several standard cavity-QED exper¬ 
iments, the atom-field coupling A is normally orders of 
magnitude smaller than the bare transition frequencies 
cj, and this makes the RWA a valid approximation 
to simplify the quantum Rabi model [3]. The straight¬ 
forward solvable nature of the Jaynes-Cummings model 
results from the presence of a continuous U{1) symme¬ 
try; the total number of excitations N = a)a is 

preserved. This should be contrasted with the discrete 
Z 2 parity symmetry of the Rabi model characterized by 
the unitary Uz 2 = exp [i (^d^a + d+d”) tt] . The action 
of Uz 2 is to flip the signs of d, d^, while leaving dz 
invariant. Thus, the Rabi model has a lower symmetry 
than the Jaynes-Cummings model, which naturally is the 
reason why finding an exact Rabi solution is very hard. 

The RWA simplified quantum Rabi model is undoubt¬ 
edly one of the most celebrated models in quantum optics 
and has received well deserved attention [30H32]. How¬ 
ever, there has been a recent surge of interest in exploring 
hybrid quantum systems in the so-called “ultra-strong” 
coupling regime [33] . In this ultra-strong coupling regime 
the atom-field coupling strength A can approach a non- 
negligible fraction of the bare transition frequencies cj, 
thereby making the Jaynes-Cummings model a non-valid 
approximation for the quantum Rabi model |34|. From 
the discussion of the introduction, it should be clear 
that in order to generate large steady state entanglement 
the coupling A should be made as large as possible (al¬ 
lowed by other appro xima tions like the two-level approx¬ 
imation). As A > y/ujU the ground state of the Rabi 
model undergoes a qualitative change [32] [35] where the 
field builds up a large non-zero photon number. In the 
Dicke model this marks the normal-sup err adiant phase 
transition [36]. This extremely large A regime has been 
termed the “deep strong” coupling regime [37]. Thus, 
as the ultra-strong coupling regime has recently been 
realized in certain (non-driven) circuit QED architec¬ 
tures [ 33 ] (A ^ cj/lO), we would like to explore the quan¬ 
tum Rabi model further, into the deep strong coupling 
regime where the coupling parameter A ^ cj, O. It is 
still unclear whether or not a non-driven system could 
attain such couplings [38]. Equally important, in this 
deep strong coupling regime the no-go theorem tells us 
that the ‘self-energy’ of the field cannot be neglected and 
by including such a term the passing to a large photon 
populated ground state will not occur [41]. Therefore, 
to overcome such hindrances we choose to work with an 
effective realization of the quantum Rabi model with suit¬ 
ably engineered Raman driving [igiiHlIlSl- Furthermore, 
as such driven models may derive more generalized ver¬ 
sions of Hrm, we will work with the anisotropic Rabi 
model (ARM) [42] 


-^ARM = wa^a+^(T2+Ai(a^(T - 
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- G d). 

(3) 


The possibility to realize the above ARM and also to 
tune co’,^J,Ai,A 2 relative to each other derives from the 
fact that amplitude of the Raman drive effectively de¬ 
termines the coupling parameters [T 5 I HHI [29] . With this 
model it is also easy to extract the importance of the 
counter rotating terms neglected in the RWA. In par¬ 
ticular, the above Hamiltonian reduces to the Jaynes- 
Cummings model when A 2 = 0, and to the quantum 
Rabi-model when Ai = A 2 = A. One important observa¬ 
tion is that the Z 2 parity symmetry is preserved for the 
anisotropic Rabi model which naturally implies that the 
eigenstates can be assigned an even or odd parity, i.e. 
simultaneous eigenstates of IJZ 2 with ±1 eigenvalues. In 
passing we may note that recently, the anisotropic Rabi 
model was also termed the [/(1)/Z2 Rabi model as it in¬ 
cludes both the Jaynes-Cummings and Rabi models as 
limiting cases [43]. A general state (not an eigenstate) 
with even(+)/odd(—) parity takes the respective forms 

|^+) = ^(4|2n + l)|t)+ct|2n)U)), 

n 

(4) 

|^_) = 5](4|2n)|t)+4|2n + l)U)), 

n 

where the first ket state |n) represents the boson Fock 
state and the second |t)/|i) fhe “up/down” eigenstate 
of the Pauli 2 ; matrix. Naturally, the coefficients cl;-^ 
and will depend on the system parameters and, due 
to the simplicity of the Rabi model, they can in princi¬ 
ple be accurately obtained by numerical diagonalization 
of (§. Let us give one more comment related to the 
symmetries of Harm- To make this comment clear we 
assume that the boson mode is only weakly populated 
such that we can truncate it to contain only the vac¬ 
uum and the first Fock state. In this case we relabel 
the creation/annihilation operators with t+/t“ respec¬ 
tively. In this truncated space we can define effective 
Pauli operators = (fa, + if^)/2, t~ = {r^ — iTy)l2^ 
and Tz = 2f+f“ — 1. Within this approximation and 
formulation we rewrite the Hamiltonian as (up to a con¬ 
stant) 


- ^ ^ Ai + A2 ^ ^ Ai 

Harm = - - - TxGx-\--\ - 


A 2 


lyUy 


• (5) 


Expressed in this form we can identify the Jaynes- 
Cummings model (A 2 = 0) with the two site Heisenberg 
XX model and the anti-Jaynes-Cummings model [44] 
(Ai = 0) also with the XX model but with one cou¬ 
pling supporting ferro- and the other anti-ferromagnetic 
order. Similarly, the Rabi model (Ai = A 2 ) can be identi¬ 
fied with the two site transverse Ising model, and finally 
the ARM (Ai 7 ^ A 2 ) with the two site Heisenberg XY 
model. Naturally, the character of the couplings are not 
just limited to the above approximation where the boson 
Hilbert space has been truncated to two states, but holds 
also for the original Hamiltonian. We cannot talk about 
phase transitions in these finite size models (contrary to 
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the Dicke model), but nevertheless the types of coupling 
share great similarities with these spin systems which all 
are critical in the thermodynamic limit [45]. We note 
that the Dicke phase transition is within the Ising model 
universality class. Even though we cannot properly take 
a thermodynamic limit of the Rabi model, it was shown 
that typical features (non-analyticity and vanishing or¬ 
der parameter in one phase) of a phase transition can 
also be achieved in the Rabi model by letting VL/uj oo. 
Not surprisingly, here one also finds universal Ising be¬ 
havior [46] . 

In m an effective optical realization of the Hamilto¬ 
nian based on multi-level atoms and cavity-mediated 
Raman transitions has been presented. Most recently 
this scheme was also experimentally demonstrated, and 
importantly the deep strong coupling regime of the Dicke 
model was reached m- A many-body generalization of 
the Hamiltonian ^ has also been recently considered, to 
explore the non-equilibrium phase diagram of the dissi¬ 
pative Rabi-Hubbard model [48|. We will not enter into 
the details of a physical realization of the Hamiltonian 
§, but will refer the reader to a comprehensive analysis 
performed in [T5] . 

Following m, damping of the field mode is the dom¬ 
inant source of dissipation and so we will approximate 
by taking it as the only dissipation channel in the open 
version of the ARM This assumption is particularly 
valid for the scheme outlined in m where the metastable 
low lying energy doublet of a Raman driven four level 
atom can act as a qubit in the ARM ® . Under the Born- 
Mar kov and secular approximations ^ | 6 ] , we model the 
evolution of the joint atom-field density operator p by 
the master equation 


P = -i[H ARM, p] + jCaP, 


(6) 


where 7 is the damping rate of the field and Cxp = 
2xpx^ — x^xp — px^x is the Lindblad super operator. It 
is worth noting that in writing the above master equa¬ 
tion we have assumed that the cavity damping is inde¬ 
pendent of the inter-mode coupling strengths Ai and A 2 . 
The form of the master equation can only arise out 
of an original time-dependent (driven) Hamiltonian. In 
the approach outlined in m the coupling Hamiltonian 
(§ is written in the frame of an external drive, imparting 
it a non-equilibrium character. For a time-independent 
system coupled to an external bath, it is required that 
the steady state will obey the principles of equilibrium 
statistical mechanics [49]. No such restriction is required 
on the dynamics arising out of a implicit time-dependent 
Hamiltonian. Naturally, this is crucial in our study, as 
we reach steady states that are different from a thermal 
equilibrium state. We thereby can conclude that even 
if deep strong coupling regime could be reached without 
external pumping such a model would be conceptually 
different from the present one as the Lindblad jump op¬ 
erators would be different [50] . 

In order to analyze the particular structure of the 
steady states pss we numerically solve the master equa¬ 
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tion and explore the resulting NESS. The steady state 
is found by integrating the equation for long times, for 
various initial states, and checking for convergence of the 
solution. According to the above discussion, a steady 
state cannot be a dark state in this case since we know 
that the dark state is the vacuum which in return is not 
an eigenstate of the Hamiltonian. We further know that 
since the system is driven, the NESS should be different 
from the system ground state which becomes the steady 
state when the system is coupled to a zero temperature 
bath [50]. However, less clear is whether the steady state 
is unique. This could, in principle, be checked by diag¬ 
onalizing the master equation, but here we only remark 
that our numerical simulations, for given parameters and 
for different initial states, all relax to steady states of a 
specific form. 

Through numerical evaluation of the master equation 
(© we arrive at the atom-field NESS density matrix pss 
and find it to have an approximate structure when order¬ 
ing the elements in a specific way. In a three-excitation 
manifold, for instance, pss can be explicitly expressed as 


Pss - 


(7) 

where the crosses mark the non-vanishing elements. The 
above structure tells us that the overlap between the 
states of opposite parities in the steady state density ma¬ 
trix Pss is identically zero. 

We can give a qualitative argument which allows the 
above structure of the density matrix pss- The Lind¬ 
blad super operator is^ obviously invariant under the ac¬ 
tion of UZ 2 ] ^aP = aU~^P' discusscd in detail 

above, the ARM © is also invariant under the action of 
parity operator Uz 2 ' Combining these two observations 
might tempt us to believe that the entire Lindblad mas¬ 
ter equation ^ is also symmetric under 1/^2 • However, 
it is important to recognize that for open quantum sys¬ 
tems a symmetry does not necessarily imply a conserved 
quantity m- Thus, in general we do not have such a 
thing as Noether’s theorem for open quantum systems. 
The numerically obtained steady state density matrix pss 
confirms this conjecture. In other words, the parity con¬ 
servation of the original ARM is broken under the 
evolution described by the Lindblad master equation ©■ 

The steady state density matrix pss 0 can also be 
expressed in a closed form as 

Pss = cos 2 6»|^'+)(^'+| +sin2 6»|^'_)(«'_|, (8) 

with the parity states Q, and 9 is some constant deter- 
mined by the system parameters. Thus, any coherence 
between the states of opposite parities is lost and the 
steady state is an incoherent mixture of these. Neverthe¬ 
less, in the steady state quantum coherence can survive 
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among the states with definite parity. The parameter 0 
is in general not 0 nor tt, meaning that parity is not a 
conserved quantity under evolution of (§); for example, 
an initial state with a definite parity will typically also 
end up in an incoherent mixture of the two parity states 
|T±). In the language of ‘pointer states’ [52], the states 
robust to the present decoherence are those with certain 
parity. 

In support of the structure of the steady state den¬ 
sity matrix (§ we argue that it is due to a symmetry 
of the equation of motion p Uz^pUz^. Let us 
assume that the steady state of the master equation (§ 
has off-diagonal coherence present between the sectors of 
opposite parities and has a structure 

/5ss = cos 2 6'|^'+)(^'+| +sin2 6»|^'_)(«'_| 

+5|«'+)(«'_|+5*|^'_)(«'+|, (9) 

where (5, likewise 6>, is some constant determined by the 
system parameters. If the above pss is a steady state of 
the master equation then it should respect the above 
symmetry i.e. 

Uz 2 PssU^^ = Pss- 

Since we get 

and we recover the above structure of the steady state 
density matrix It is worth pointing out that a two- 
site steady state density matrix of the dissipative trans¬ 
verse field Ising model can be deduced as a special case 
of our density matrix structures 0,® m- This is not 
a surprise since both the Rabi and the transverse field 
Ising models have discrete symmetries and which are also 
obeyed by the respective equations of motion. 

As seen from from Eq. the presence of the driv¬ 
ing together with the reservoir demolishes any quantum 
coherence between the two different parity components. 
Nevertheless, the steady state will in general be an in¬ 
separable state of the atom and the field. To quantify 
the amount of this bi-partite entanglement present in the 
atom-field joint density matrix we compute the logarith¬ 
mic negativity defined as log(5]]^ where l^il is 
the sum of the absolute values of all the eigenvalues of 
the partially transposed density matrix [53j |54] . A non¬ 
zero value of the logarithmic negativity is sufficient to 
ensure inseparability of the steady state density matrix 
Pss- The steady state logarithmic negativity, obtained 
from numerical integration of the master equation to first 
obtain pss and then partially transpose it, is plotted as 
a function of the dimensionless coupling ratio A 2 /A 1 in 
Fig.B As can be seen, the steady state entanglement 
between the atom and the field grows with the value of 
A 2 /A 1 . In the limit A 2 /A 1 = 0 the case of the Jaynes- 
Cummings model is recovered and we know that the 
steady state is simply p^s = |0)(0| (8) | |. Thus, the 

counter rotating terms are responsible for the build-up 


of photonic and atomic excitations and thereby also for 
bi-partite entanglement in the system. However, in the 
opposite limit, A 1 /A 2 = 0, we obtain the anti-Jaynes- 
Cummings model and here the steady state is again sep¬ 
arable; pss = |0)(0| (8) I t)(t I- A maximum of the entan¬ 
glement should therefore be given for some finite A 2 /A 1 , 
and from Fig.j^we see that it happens at A 2 /A 1 ~ 1.6, i.e. 
it is favourable to support a stronger coupling A 2 in order 
to achieve a non-classical entangled state. Even though 
the Jaynes-Cummings and anti-Jaynes-Cummings mod¬ 
els are mathematically equivalent (in the bare basis, for 
example, the Hamiltonian maintains a block structure 
with 2x2 blocks apart from the ground state), we see 
that the presence of the Lindblad term jCa alters this sym¬ 
metry, i.e. the maximum entanglement is not obtained 
for Ai = A 2 . 



FIG. 1: (Color online) Steady state entanglement (logarith¬ 
mic negativity) in the joint atom-field steady state pss, plotted 
here as a function of the dimensionless coupling ratio A 2 /A 1 . 
In the limit of vanishing Ai or A 2 the steady state is sepa¬ 
rable. The maximum entanglement obtained for A 2 ~ I. 6 A 1 
demonstrates that the presence of photon losses break the 
symmetry between the Jaynes-Cummings and anti-Jaynes- 
Cummings models. In particular, the counter rotating terms 
are responsible for counteracting the losses. The other dimen¬ 
sionless parameters are Q = a; = 1 and 7 = 0 . 1 , and for the 
plot we fix Ai = 0.5. 

The reduced density matrix of the boson field is ob¬ 
tained by tracing over the atomic degrees of freedom, 
/Sfieid = Tratom [Pss], giving 

CX3 CX3 

Pfield ~ E E P(n, m)|n)(m| + h.c. . (10) 

n=0 m=n,n+2... 

As a result of the mixture 0, the steady state field dis¬ 
tribution is an individual coherent mixture of even and 
odd numbered Fock states. However, these two coherent 
mixtures exist in two different ‘sectors’ with no coherent 
overlap between them. In Fig. we plot the steady state 
probability distribution P{n) = (n|pfieid|^) for three dif¬ 
ferent values of the ratio A 2 /A 1 , and again the steady 
state has been found by integration of Eq. 0 - We find 
that the steady state density matrix pfieid follows super- 
Poissonian statistics where the difference between the 
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FIG. 2: (Color online) Steady state photon probability dis¬ 
tribution P(n) for different coupling strengths Ai and A 2 . It 
is noted that the distribution is found to be super-Poissonian 
in all cases, i.e. the photon distributions alone does not dis¬ 
play non-classical features. The lines between the points are 
for guiding the eye, and the dimensionless parameters are the 
same as in Fig. ^ Q = ct; = 1 and 7 = 0.1. 


variance and the mean of the photon number distribution 
increases with the value of the ratio A 2 /A 1 . In particu¬ 
lar, and as expected, for increasing atom-field coupling 
the mean number of photons increases. The counter ro¬ 
tating terms play a more important role for this increase 
to occur. To explore the importance of the bare state 
coherences in Fig. we plot the absolute value of the 
first four leading order off-diagonal components of pfieid, 
namely |P(n,m)| = |(n|pfieid|^)|- It is clear that these 
imprint the coherence among the bare states belonging 
to either even or odd parity sectors. 
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FIG. 3: (Color online) Leading off-diagonal elements of the 
reduced density matrix pfieid ( |10| ), when A2/A1 = 0.5 (top 
row), A2/A1 = 1.0 (middle rov^ and A2/A1 = 2.0 (bottom 
row), and Ai = 0.5,0.75,1.0 (increasing from left to right in 
all three rows). The alternating zero and non-zero elements 
reflects the population of states with different parities. The 
other parameters are the same as for Fig. 


So far we have seen that for the driven/dissipative 
Rabi model non-classical features may survive. An in¬ 
teresting question then raises, which of the features of 
the (equilibrium) ground state of the Rabi model ^ can 
be reproduced in the ^on-equilibrium) steady state of 
the master equation (1^? For the Rabi model, in the 
deep strong coupling regime the ground state becomes a 


Schrodinger cat [32l [35] . The more resolved the cat state 
becomes (the smaller the overlap between the two sepa¬ 
rated field states in phase space becomes), the higher the 
atom-field entanglement obtained. More precisely, these 
‘dead’ and ‘alive’ field states measure the dipole moment 
of the atom (corresponding to eigenstates of ax). Thus, 
when tracing out the atom, the field state will become a 
statistical mixture between ‘dead’ and ‘alive’. Note that 
pure cat state of the form |CAT±) oc (|a) ± | — a)), for 
coherent states | ±a) with \a\ >0, have a definite parity 
and the corresponding photon distributions contain ei¬ 
ther only even or odd photon states |n). Thus, we ask if 
we can expect a similar statistical cat also for the reduced 
density matrix for the boson field ([1Q|)? Figure [^reveals 
at least that the photon distribution of the steady state 
is super-Poissonian, and combining this with Fig. one 
may expect the cat structure to survive photon decay. 
This should be contrasted with the pioneering experi¬ 
ments in the ENS Paris group which measured the decay 
of the cat into a fully separable atom-field state before 
all photons had leaked the cavity m- The phase space 
distribution of the steady state field distribution for the 
Rabi model is indeed split in two as is shown in Fig. 
The Wigner function for the field mode can be obtained 
from the following transformation of the density operator 


W(a) = -TT[D\a)pfi,iaD{a){-lf% (11) 

TT 


where D{a) = exp((adl — (a*d) is the displacement opera¬ 
tor and a = ar E ioLm is a complex parameter |56]. Note 
that for other expressions for the Wigner function, the 
real and imaginary parts of a are related to ‘momentum’ 
p and ‘position’ x. As can be seen from the figure, on 
increasing the value of the ratio A 2 /A 1 the separation be¬ 
tween the two-lobes becomes more pronounced and even¬ 
tually the “two-lobes” cease to overlap with each other. 
A similar two-lobe structure has been found in the tran¬ 
sient dynamics of the Wigner function for the dissipative 
quantum Rabi model with an additional Stark shift term 
included [29]. For the closed system, this behavior readily 
follows from a sort of mean-field approach identical to the 
Born-Oppenheimer approximation [32l|57]. In particular, 
within this approximation and for 7 = 0 the maxima of 
the Wigner function are found for 


Im((a) 


0 , 



9 < 9c 

( 12 ) 

9 > 9c 


where g = (Ai + A 2)/2 and Qc = \/<^r^/2, and Re(a) = 0. 
Note that Qc coincide with the critical coupling of the 
Dicke model [36]. Thus, the splitting grows approxi¬ 
mately linear with Ai + A 2 . If the signs of the couplings 
Ai and A 2 are different the roles of the real and imag¬ 
inary parts of a are interchanged. Now, Fig. [^ shows 
that the real parts for the locations of the Wigner func¬ 
tion maxima are actually non-zero. This is seen in the 
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anti-clockwise tilting of the Wigner function. This is a 
result of the Lamb shift [56], i.e. it only arises due to 
the coupling to the bath. This is also easily understood 
by considering the mean-field versions of the Heisenberg 
equations of motion for the open Rabi model. The steady 
state solution for the field will then include a non-zero 
imaginary part (or real part depending on how one de¬ 
fines a) which implies that the field phase for the sepa¬ 
rated blobs is not exactly Ttt [58] . 



Re(o) Re(o) Re(o) 

FIG. 4: (Color online) Steady state Wigner function 0 for 
the anisotropic Rabi model; A 2 /A 1 = 0.5 (top row), Rabi 
model A2/A1 = 1.0 (middle row), and A2/A1 = 2.0 (bottom 
row), and Ai = 0.5,0.75,1.0 (increasing from left to right in 
all three rows). The splitting between the two blobs scale as 
Ai + A2 for the closed case. The presence of photon loses tend 
to make the splitting smaller since the corresponding Lindblad 
term favor a vacuum state. The bath also brings about an 
energy (Lamb) shift which is manifested in the tilting of the 
Wigner function. The absence of negativity of the Wigner 
functions in this plot does not solely result from the reservoir 
induced decoherence but also from the entanglement shared 
between the atom and the field. The other parameters are 
the same as for Fig. 


III. MODEL B - ANISOTROPIC DICKE MODEL 
(NORMAL PHASE) 

The Rabi model is highly quantum in the sense that 
quantum fluctuations play a dominant role in the spin 
sub-space, and consequently it lacks a natural classical 
limit. Such a limit would, however, emerge if we allow 
the spin to become very large, i.e. the Pauli matrices 
are replaced by general angular momentum operators Sa 
{a = X, z). In doing this replacement we obtain the 
anisotropic Dicke model (ADM) 

HADM=toa^a^S,+Xi S-+S+aj+X2(a^S++3-0^. 

(13) 

It is the purpose of this section to demonstrate that non- 
classicality of the steady state survives also in this clas¬ 
sical limit. 

With the correct scaling of the coupling parameters 
(relative to the system size), the Dicke model is critical 


with the normal phase characterized by a vacuum photon 
field and the spin pointing towards the south pole. The 
superradiant phase instead comprises a coherent photon 
state and the spin rotated away from the south pole [36] . 
True criticality only emerges in the thermodynamic limit 
meaning that the spin S ^ oo ED]. In the strict limit 
of infinite spin, the spectrum is linear and we may re¬ 
place spin operators by boson operators (plus an overall 
energy shift). In such a case we consider the following 
Hamiltonian 


^nAD = ood^d QPb -f Ai (^d^b + Pd^ + A2 + bd^ , 

(14) 


where the notation ^nAD refers to the linearized 
anisotropic Dicke model in the normal phase. Clearly 
the above Hamiltonian can be used to model various 
other physical systems, including coupled harmonic os¬ 
cillators/cavity modes [49]. Going from the non-linear 
model (13) to the linear approximated model (14) al¬ 
lows for an analytical approach. For finite systems the 
non-linearity will enter after sufficiently long times [6T] . 
Here, however, we assume that this time-scale is much 
longer than the intrinsic relaxation time of the open sys¬ 
tem, i.e. we push any quantum revivals to happen be¬ 
yond our experimental times. In the normal phase of 
the Dicke model, the linear Hamiltonian HuAB describes 
the collective quantum excitations/fluctuations (Bogoli- 
ubov modes) [15] . This is most easily seen by linearizing 
the original Dicke model in the Holstein-Primakoff boson 
representation [niEn]- In the superradiant phase, the re¬ 
sulting quadratic boson Hamiltonian contains additional 
squeezing terms, e.g. 6^ and [62]. Within the effective 
linear model, the quantum fluctuations cause entangle¬ 
ment in the system, and for the closed Dicke model it 
has been found that the atom-field entanglement peaks 
at the critical point [63]. We point out, however, that 


since the Lindblad jump operator [d, ^adm] 7 ^ 0 it is not, 
in principle, clear that criticality of p^s is preserved in the 
presence of the photon decay. This issue was recently ad¬ 
dressed and photon losses do not destroy criticality but 
alters both the critical coupling and the critical expo¬ 
nents [151 [ 59 ]. As discussed in more details below, the 
conservation of the critical point under photon losses can 
be understood from the fact that the state for the normal 
phase is actually a dark state of the Lindblad operators. 
Here, as we are interested in the opposite (classical) limit 
of the Rabi model, we only consider the fluctuations in 
the normal phase and not in the superradiant one. This 
emerges in the present analysis as destabilization of the 
Bogoliubov modes upon approaching the critical point. 

Going beyond the study in [63], we here include the 
interaction with the outside environment. We again en¬ 
visage a scenario where damping of the field mode is the 
only dominant dissipation channel; 


P = -i[HnAD, p] + l^aP- 


(15) 


In the above master equation we have neglected dissipa¬ 
tion of the mode b. We again argue that this assump- 
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FIG. 5: Stability criterion ( 0 for different values of the 
ratio A 2 /A 1 , shown here as a function of the decay rate of 
the field mode 7 and the coupling strength Ai . The non-zero 
value C 7 ^ 0 mark the transition to the unstable regime. The 
dimensionless atomic and field frequencies Q = cc = 1 . 


tion can be justified if it is possible to realize the ADM 
(14) through a scheme outlined in m- Specifically, mode 
b can play the role of collective Bogoliubov excitations 
of a collection of “two-level” systems, where each “two- 
level” system can be the metastable low lying energy dou¬ 
blet of a Raman driven four level atom. Also, as discussed 
in the previous section, a phenomenological introduction 
of the Lindblad super operator for the field mode as in¬ 
troduced in the master equation (15) can only be jus¬ 
tified for an original time-dependent Hamiltonian. It is 


worth pointing out that for a time-independent model 
described by the Hamiltonian (14) a correct description 
of the open dynamics results in a non-local master equa¬ 
tion which disagrees with the predictions of the master 
equation (15) [49] . 

As already mentioned, since the Hamiltonian is 
quadratic we may approach the problem analytically. To 
this end we define the two-mode normal-ordered char¬ 


acteristic function x(ea,e 5 ,t) = e 


ebb^ ell 


')• 


From the quantum characteristic function the complete 
statistical description of the corresponding state can be 
obtained. One can therefore obtain the expectation val¬ 
ues of quantum mechanical observables, e.g. for the equal 
time correlator one has 

= (^) (^) 

For an initial Gaussian state of the two modes d and 6 , 
the master equation (15) can be expressed as a partial 


differential equation for % 


d 


FIG. 6 : Bi-partite entanglement between the modes a and 
b for three different values A 2 /A 1 and as a function of 7 and 
Ai. Right before the solution gets unstable (white region, and 
see previous figure]^ the entanglement is maximized. This is 
in agreement with the fact that in the full Dicke model this 
marks a critical point m- Note, however, that the critical 
point has been shifted due to the coupling to the bath. The 
remaining dimensionless parameters are Q = ct; = 1 . 


where z'^ = {ea, el,eb, ej), V = (gf^, and 

iA2/2 

n n 

M = 


N = 


are the drift and diffusion matrices respectively. We first 
analyze the steady state solution of the master equa¬ 
tion (15). In particular, the solution is stable if the real 
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parts of all the eigenvalues of the diffusion matrix N are 
all negative. In particular, when = cj, the stability of 
the steady state requires 


( = Max ±iy ) — 7 < 0, 


(17) 


where. 


/+ — 


■ 7 ^ 0 ;^ 


= z'^Mzx + 2^NVx, 


(16) 


= ^- 4(^X1-XI+u:^ + ^ 4Xloj^ 

v~ = ^ 7^—4 — X\ 4- — \J4X^10^ — . 

The stability parameter C is shown in Fig. as a func¬ 
tion of the damping rate of the field mode 7 and the two¬ 
mode coupling strength Ai. As is evident from the figure, 
increasing the damping of the field mode allows a stable 
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steady state to be achieved for an extended region of pa¬ 
rameter space spanned by Ai. This behavior can be un¬ 
derstood physically by considering the different terms of 
the master equation (15). First we may note that the in¬ 


stability of the mode coincides with the Dicke phase tran¬ 
sition. Indeed, the model Hamiltonian (14) describes the 


collective excitations in the normal phase and not in the 
superradiant phase, and hence, these excitation modes 
cannot be analytically continued into the superradiant 
phase. Thus, we conclude that increasing the photon loss 
gives a more extended normal phase. This is expected 
since the Lindblad term of Eq. (15) favors the vacuum 


state of the field, which is also favored by the bare field 
energy ujd)a. On the other hand, the atom-field coupling 
term is not minimized by such a field state. Without pho¬ 
ton decay, the transition occurs at (Ai + A 2 ) = \/ujVt (see 
Eq. (p!^), but now losses support the normal phase and 
thereby shift the critical value. Eollowing Refs. p!5l [64) 
one explicitly finds Ac = + 7 ^) /2 which is ob¬ 

tained at the mean-field level, i.e. in the thermodynamic 
limit. In the Dicke limit (Ai = A 2 = A) it is easy to obtain 
the phase boundary + 7 ^ = 4A^, which retrieves the 
result for the critical coupling of the open Dicke model 
on resonance [151 El • 

The steady state of light-matter coupled quantum sys¬ 
tem modeled under the master equation was shown 
to be an inseparable state of the cavity field and the 
two-level system. Somewhat surprisingly, this observa¬ 
tion can also be extended to the steady state of the 
normal phase of the Dicke model (15). In particu¬ 


lar, one could argue that in the thermodynamic limit 
as studied here any entanglement should vanish since 
for the Dicke model quantum fluctuations are negligi¬ 
ble [65] . That this is not the case derives from the pres¬ 
ence of the critical point. To demonstrate sustainable 
quantum correlations we next evaluate the steady state 
bi-partite entanglement between the two modes a and 
b and characterize it in terms of the logarithmic neg¬ 
ativity which for two-mode Gaussian states serves as 
a necessary and sufficient criterion for the inseparabil- 
ity [ 66 ] . A two-mode Gaussian state can be fully quanti¬ 
fied in terms of its covariance matrix V which is a 4 x 4 
symmetric matrix with Vij = {{RiRj + RjRi))/2 and 
= {Qa,Pa,qb,Pb)- Here qa,b and Pa,b are the posi- 
tion and momentum quadratures of the mode a{b). Eor 
a two-mode Gaussian continuous-variable state with co- 
variance matrix V, the logarithmic negativity is obtained 
as A/' = Max[0, —log( 2 z^_)] [ 66 ], where U- is the smallest 
of the symplectic eigenvalues of the covariance matrix, 

given by U- = a/2 — (a^ — 4DetV)/2. Here 

a = DetAi + DetBi — 2DetCi 


V = 


Ai Cl 


where Ai (Bi) accounts for the local variances of mode 
a (b) and Ci for the inter-mode correlations. Using the 
numerical solutions of the partial differential equations 


(16) we compute the logarithmic negativity. The steady 
state two-mode entanglement is shown in Eig. [^ where 
it is given as a function of the damping rate 7 and the 
two-mode coupling strength Ai. As is evident, the maxi¬ 
mum attainable bi-partite entanglement between the two 
modes monotonically increases on increasing the value 
of the ratio A 2 /A 1 . The highest entanglement is ob¬ 
tained at the critical point, which is also known from 
the closed Dicke model [63]. This is, in fact, a general 
property for quantum phase transitions; at the critical 
point where characteristic length scales diverge, the en¬ 
tanglement also diverges m- However, already for the 
closed Dicke model the phase transition is not a ‘typical’ 
quantum phase transition since there is no length scale 
in the problem and, as we mentioned above, quantum 
fluctuations vanish in the thermodynamic limit. Because 
of this, the transition has been called ‘classical’ [65]. As 
such, the results of Eig. [^suggest that, qualitatively, the 
same behavior of the entanglement is found for a dynam¬ 
ical phase transition of an open driven system. 


IV. QUANTUM CONTROL IN HYBRID 
ARCHITECTURES 


i __ ._i 


Source Cavity 


Driven Cavity 



i 



FIG. 7: Schematic plot of the feedback loop between the two 
cavities. We envision a scenario of two coupled cavities; a 
source cavity and a driven cavity. The driven cavity is as¬ 
sumed to be slaved to the source cavity. The source cavity 
drives the state of the slave cavity by means of a unidirec¬ 
tional coupling (shown by a solid line). The driven cavity in 
turn also influences the state of the source cavity by means 
of a reversible interaction (shown by a dotted line). The si¬ 
multaneous presence of reversible and irreversible couplings 
between the cavity modes results in an all-optical feedback 
loop. 

In the previous two sections we have studied the open 
dynamics of the anisotropic Rabi and Dicke models and 
have explicitly explored the properties of the NESSs of 
these models. However, so far we have eluded ourselves 
from effectively engineering the open dynamics of these 
hybrid models. Even though it was found that entangle¬ 
ment persists in both models despite losses, in the present 
section we wish to explore the possibilities to enhance the 
amount of entanglement by actively introducing feedback 
into the model. 

An unusual ingredient which enriches the dynamics of 
quantum systems is their inherent ‘openness’ which may 
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FIG. 8: Illustration of the general all-optical feedback of Fig.|^ 
specifically used for establishing quantum control in the Dicke 
model of Section [ni| T he internal dynamics of the source cav¬ 
ity is modeled by (14) while the driven cavity is initially pre¬ 
pared in its vacuum state. Mode a of the source cavity and 
the mode c of the driven cavity are interacting under a re¬ 
versible interaction of the form (18). Mode a of the source 
cavity also couples irreversibly to mode c of the driven cav¬ 
ity. Such a feed-forward coupling can be engineered through 
feeding the mode d from one end of the source cavity and 
reflecting it onto the driven cavity through a series of mirrors 
(filled) and beam-splitters (unfilled). Optical circulators or 
Faraday isolators can be used to prevent interference from re¬ 
flections in the opposite directions. Reversible and irreversible 
interactions jointly constitute an all-optical feedback. If this 
feedback loop has negligible time delay a Markovian master 
equation to describe the open dynamics of the source cavity 


can be derived (21) as detailed in Section IV A 


result in exotic NESS properties [3 [m [48l |68l [69] . This 
was in particular the theme of the previous two sections. 
However, engineering a controlled degree of dissipation 
is also an important tool in realizing tasks of quantum 
information processing. In this spirit we will now use an 
all-optical feedback scheme for quantum state protection 
and establishing quantum control in the normal phase 
of the Dicke model. The scheme we choose is based on 
coherent feedback, i.e. no measurement is performed in 
the feedback loop m- A symbolic plot of a coherent 
(all-optical) feedback scheme is illustrated in Fig. it 
comprises of two units, which have been labeled source 
and driven cavities, with reversible and irreversible cou¬ 
plings between them. The simplest all optical feedback 
loop will use the output from a source cavity and feedfor¬ 
ward into a driven cavity which is coupled to the source 
cavity in some way. It is worth mentioning that unidirec¬ 
tional and bidirectional couplings between the source and 
the driven cavities jointly constitute an all optical feed¬ 
back scheme. A coherent feedback loop, however, can 
also introduce a finite time-delay r. Firstly, assuming 
that the time-delay r introduced by the feedback loop is 
negligible, we study the open dynamics of the source cav¬ 
ity P^H2Q] . Subsequently, making use of a time-delayed 
feedback control method, [211423] we incorporate a time- 
delay r introduced by the feedback loop. 


A. All optical feedback loop with negligible 
time-delay r 


We now apply the general scheme of all-optical feed¬ 
back illustrated in Fig. for a specific task of quantum 
state protection in the anisotropic Dicke model of Sec. 


Ill 


__ A schematic of our proposal to implement all-optical 

feedback is shown in Fig. The internal (unitary) dy¬ 
namics of the source cavity is described by the Hamilto¬ 
nian (14), i.e. ^source = ^nAD and the driven cavity is 


reversibly coupled to the field mode of the source cavity 
by the following Hamiltonian UHlllg] 




IjJ, 




(at 


ta), 


(18) 


where /x is a dimensionless coupling parameter, ct and c 
are the creation and annihilation operators for the field 
mode of the driven cavity which has damping rate 7d. 
The driven cavity is assumed to be prepared in its vac- 
uuni state with its respective internal dynamics modeled 
as ^driven = ^d^^c. A revcrsiblc interaction of the form 


(18) can arise through mode overlap between the modes 
d and c [70] . 

Under the Born-Markov approximation a joint state of 
the source and driven cavities, represented here as FF, 
evolves under the following master equation [SI nil US] 


w = -i 


47nAD T 47driven T 


+V^{[aW,c^ + [c,Wa^) 

+ '^CiiW +'^CeW. 

2 2 


(19) 


It should be remarked that W represents a tri-partite 
state of modes d, 6, and c. On the grounds of arguments 
presented previously, we have again neglected damping 
of the mode b of the source cavity. The two terms ap¬ 
pearing in the second line account for the unidirectional 
coupling between the source and driven cavities and the 
last two terms are the individual Lindblad operators de¬ 
scribing photon losses for the source and driven cavity 
modes respectively. As illustrated in Fig. such a uni¬ 
directional coupling between the source and driven cavi¬ 
ties can be established using an optical circulator [18]: a 
non-reciprocal optical device such as Faraday rotator can 
be used to establish irreversible coupling between the op¬ 
tical modes d and c of the source and the driven cavities. 
As pointed out in Refs. miiis], for the feedback loop to 
be effective the driven cavity should respond much faster 
than the source cavity. We thus work in a regime where 
7d ^ 1 meaning that the state of the driven cavity is 
slaved to the state of the source cavity. Since 7 ^ is the 
large parameter (i.e. determining the fast time scale) and 
the driven cavity is assumed to couple to a zero temper¬ 
ature reservoir, the population of the driven cavity mode 
can be truncated to only the lowest photon states. We 
therefore approximate the joint state of the two coupled 
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cavities as 




Poo|0)(0| + /3io|l)(0| + /3|o|0)(l| 

+Pii|l)(l| + P2 o|2)(0| + pIq\0){2\, 


( 20 ) 


where pij is the conditional state of the source cavity 
(a joint state of modes a and b) when the driven cavity 
is projected on the state space \i){j\. Using the above 
ansatz and adiabatically eliminating the driven cavity it 
is possible to derive an effective master equation for the 
source cavity alone [TSl |T9] . Following the derivation pro¬ 
vided in Appendix one arrives at the resulting master 
equation 


P = Poo + Pii = — 


HnAB.P 


7eff p ^ 


( 21 ) 


where 7 eff = 7 (1 H- /i (2 -|- /i)) is the new effective damp¬ 
ing rate of the field mode. On choosing the dimensionless 
parameter p = — 1 , it is remarkable to observe that an 
all optical feedback loop is capable of completely block¬ 
ing the loss of the source cavity. However, losses in 
the feedback loop would deteriorate the effectiveness of 
such feedback. If the feedback loop has an efficiency 
7 (< 1)5 then it is possible to show that the effec¬ 
tive damping rate of the source cavity gets modified as 
7eff =7(1+W(2 + Ai)) US]. We therefore conclude that, 
under the assumptions that the time-delay introduced by 
a feedback loop is negligible and the state of the driven 
cavity is effectively slaved to the source cavity, it is pos¬ 
sible to establish arbitrary control over the damping rate 
of the source cavity. This can be an important step for 
quantum state protection in hybrid quantum systems. 

As an example application of an all optical feedback 
scheme for quantum state protection, we assume that the 
two coupled modes interacting under the Hamiltonian 
(14) are initially prepared in a NOON state [71]. NOON 
states are Bell-like entangled states with applications in 
quantum metrology and quantum lithography and they 
may also be utilized for achieving phase supersensitivity 
[72] . We assume that the modes a and b are initialized 
in a state |T^) = (|A^)a|0)6 — \0)aN)l,)/^/2. We examine 
the autocorrelation function, i.e. the overlap between 
p(0) and the time evolving joint state of the two modes 
p(t), which for mixed states is given by the Uhlmann 
fidelity [73l [74] 


$ = Tr^ v7W/5(t) vTM- 


( 22 ) 


The time evolution of the fidelity 4> for different values 
of N and two different values of 7eff is plotted in Fig. [^ 
(a) and (b). We note that the results are for a small 
coupling ( 2 OA 1 = 2 OA 2 = uj = ft) in order to assure nu¬ 
merical convergence. In Fig.|^(c) we also show the time 
evolution of the fidelity 4> for an initial bi-model NOON 
state in absence of photon losses. As is evident, by con¬ 
trolling the damping rate of the source field it is possible 
to prolong the lifetime of the initially prepared NOON 
state. It follows that suitably choosing 7eff allows the 
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FIG. 9: (Color online) The upper plots (a) and (b) give the 
time evolution of the fidelity 4> for a bi-modal NOON state 
in the presence of losses, while (c) shows the same but in ab¬ 
sence of photon losses. By comparing the first two examples 
it is clear that the feedback loop makes it possible to preserve 
coherent evolution for longer times. The last plot (d) is the 
same but for an entangled cat state with one mode in vacuum 
and the other in a coherent state ja). The shorter period re¬ 
sults from the dephasing of the different Fock states involved. 
The other dimensionless parameters are U = u = l,Ai = 0.05 
and A 2 /A 1 = 1 in all four panels. 


value of the fidelity to be kept above its classical value 
of 2/3 [ 75 ] for a longer time. It is also worth noting that 
NOON states become more susceptible to the environ¬ 
mental damping with the increase in the excitation num¬ 
ber N. In this scenario establishing a control over the 
damping rate of the field in the source cavity could be an 
important step in preserving initial quantum coherence. 
In the Appendix we analytically analyze the structure 
of especially plot (c), and explain why the fidelity stays 
constant for A" = 1, 2 and not for A = 3, 4. 

Also shown in Fig. (d) is the evolution of the 
Uhlmann fidelity of a two-mode entangled coherent state 
1^') = (|a)a|0)b + |0)a|a)6)/y2(l + e-l“l^), where |q;) 
represents a coherent state. Entangled coherent states 
can be a useful resource in quantum metrology and can 
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exhibit noticeable improved sensitivity for phase estima¬ 
tion when compared to that for NOON states. Entangled 
coherent states can also outperform the phase enhance¬ 
ment achieved by NOON states both in the lossless, weak, 
moderate and high loss regimes US]. As demonstrated in 
Fig. [^(d), when compared to a NOON state a two-mode 
entangled coherent state is found to be more resilient to 
photon losses. Nevertheless, the quantum fidelity of en¬ 
tangled coherent states also seems to decay significantly 
with the increase in the mean number of photons 


B. All optical feedback loop with finite time-delay 


As mentioned before, in deriving the master equation 
it has been assumed that the time-delay r intro¬ 
duced by the feedback loop is negligible [18]. In this 
section we explore a complementary regime when the 
feedback loop introduces a non-negligible time-delay r. 
In particular, we use the time-delayed feedback control 
method of Ref. [2TU23] and apply it to our Hamiltonian 
( [l^ . We refer the reader to m for a specific proposal 
implementing a finite time-delay in an optical feedback 
loop and applying it to the Dicke model ( [m]). In this 
work we go beyond the analysis presented in j23j and will 
include quantum fluctuations to check the steady state 
stability of the time-delayed feedback control scheme and 
to compute steady state correlations between modes d, b. 
Considering a specific all-optical time-delayed feedback 
control strategy discussed in detail in [23] and apply¬ 
ing it to our Hamiltonian (in we arrive at the following 
Heisenberg-Langevin equations of motion (14) 


= AV{tf - 


+ BY{t-Tf, (23) 


where 

V(i) = {a{t),a\t),b{t),b\t)), 


Langevin equations we have assumed that the time- 
delayed feedback loop has unit efficiency [23]. To con¬ 
nect with the approach of the previous sections, we have 
again assumed that the damping of the mode a is the only 
dominant channel of dissipation. Also, the time-delayed 
control feedback strategy (23) is implemented through a 


control force which is generated from the difference be¬ 
tween the instantaneous cavity field a{t) and the field at 
some point in the past a{t — r) [23] . 

As a first step to check the influence of the time-delayed 
feedback control on our hybrid quantum system we check 
the stability of the above Heisenberg-Langevin equations 
semi-classically, i.e. we assume the noise (Vin(^))=0. Us¬ 
ing an ansatz V{t) ^ we get the following secular 
equation mjna 


det(A-B + Be-'^^ - Al). 


(25) 


For a non-zero value of r, this transcendental equation 
has an infinite set of complex solutions for eigenvalues A. 
The steady state is stable only if the real parts of all the 
solutions A are negative mum- We numerically solve 
the above secular equation ( [ 2 ^ for all possible roots A. 
In doing so we choose Ai and 7 such that in the absence 
of a time-delayed feedback loop the steady state is stable 
for all values of the ratio 0 < A 2 /A 1 < 2, see Fig. In 
Fig. we plot the maximum of the real part of all possi¬ 
ble solutions A. We find that at a semiclassical level, and 
for the set of parameters considered in Fig. our time- 
delayed feedback control strategy does not qualitatively 
alter the stability of the steady state. 



A = 



= ia{t - 

- T),a^t 

-t)'. 

,b{t - 

) = 

(a-(i), 

at“(i),0 

',0), 


/- 

-iuj — 7 

0 

-iAi 

—iX[ 

0 

iuj — j 

to 

iX\ 


—iXi 

—i\2 

— iU 

0 


to 

i\i 

0 

iQ 


B = 


/ 7/2 0 0 0 

0 7/200 

0 0 0 0 

Vo 0 0 0 


(24) 

Here, Vin{t) contains the input noise terms [77] • It should 
be pointed out that in writing the above Heisenberg- 


FIG. 10: (Color online) Maximum of the real part of all pos¬ 
sible solutions of the secular equation ( |25| ) for three different 
values of the ratio A2/A1 and shown here as a function of 
the time-delay r. The dimensionless parameters have been 
taken Q = u = I,Ai = O.I and 7 = O.I. For these set of 
parameters the steady state is stable for all values of the ra¬ 
tio 0 < A2/A1 <2 in the absence of a time-delayed feedback 
loop, see Fig. 

To show that our time-delayed feedback control strat¬ 
egy is indeed capable of influencing the steady state be¬ 
haviour of our hybrid quantum system, we now embark 
on a full quantum treatment to explore the steady state 
correlations between the modes a and b with their dy¬ 
namics governed by Eq. ( [ 2 ^ . We solve the equations in 
the Fourier space and reconstruct the steady state two- 
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FIG. 11: (Color online) The upper plot (a) displays the bi¬ 
partite entanglement (negativity) between the modes a and b 
as a function of the dimensionless ratio A 2 /A 1 and the time- 
delay T, and in the lower plot (b) the entanglement is given 
for different ratios A 2 /A 1 and for different choices of r. Also 
shown in (b) is the steady state bi-partite entanglement be¬ 
tween the modes a and b in the absence of a time-delayed 
feedback loop, and we especially note that by introducing a 
time-delay the amount of entanglement can be increased. In 
the black region in (a) the solution does not approach a phys¬ 
ical steady state. The dimensionless parameters have been 
taken Q = u = l,Ai =0.1 and 7 = 0.1. 


mode covariance matrix in order to capture the quan¬ 
tum correlations between the modes d and b. We then 
evaluate the bi-partite entanglement (logarithmic nega¬ 
tivity) between the two modes d and b. In Fig. [^(a) we 
display the bi-partite entanglement between the modes 
d and 6 as a function of the dimensionless ratio A 2 /A 1 
and the time-delay r. As for Fig. we choose Ai and 
7 such that in the absence of a time-delayed feedback 
loop the steady state is stable for all values of the ratio 
0 < A 2 /A 1 <2. We find that for a fixed value of the time- 
delay r the bi-partite entanglement between the modes 
d and b increases with the value of the ratio A 2 /A 1 , this 
is expected since the coupling between the two modes is 
also increased. However, we also observe that there is a 
“delay-window” (shown in black and labeled ‘Unphysi- 
caF in Fig. 11 (a)), such that if r is chosen within this 
interval then a physical steady state is never reached [79] . 
We also find that a suitable choice of the time-delay can 
also increase the steady state bi-partite entanglement be¬ 
tween the modes d and b. More precisely, the bi-partite 
entanglement between the modes d and b for different 


choices of r is given in Fig. 11 (b), where for the sake 
of elucidating the role of time-delayed feedback loop we 
have also shown the steady state entanglement between 
the modes d and b in the absence of a time-delayed feed¬ 
back loop. 


V. CONCLUSIONS 

In this work we have studied two topical hybrid sys¬ 
tems, the anisotropic Rabi and Dicke models, both of 
which combine disparate quantum degrees of freedom; 
light and matter. While similar, the two models oper¬ 
ate in different ‘regimes’; the Rabi model displays large 
quantum fluctuations which is not the case of the Dicke 
model. We have first explored the open dynamics of 
a two-level system (a model atom) strongly coupled to 
a boson (photon) field, i.e. the Rabi model. We have 
shown that the NESS is an entangled state of the light 
and atom systems. More precisely, the steady state is a 
statistical mixture of two states with opposite parities. 
In phase space these even and odd parity states forms 
a cat-like structure, which however gets ‘tilted’ due to 
the reservoir-induced Lamb shift. As a second physical 
model we have investigated two coupled boson modes 
which can describe linearized interactions between a cav¬ 
ity field and a collection of two-level systems in the ther¬ 
modynamic limit, i.e. the Dicke model. This linearized 
version of the anisotropic Dicke model describes the col¬ 
lective excitations in the normal phase. We have espe¬ 
cially explored the stability and bi-partite entanglement 
in this second physical setting, again in the presence of 
coupling to an environment. Even though quantum fluc¬ 
tuations are negligible in the thermodynamic limit, en¬ 
tanglement survives which we attribute to the presence 
of a critical point which persists also under dissipation. 
As a way to establish quantum control and increase non- 
classical properties in such hybrid quantum architectures, 
we propose to use an all optical feedback strategy. We 
demonstrated this approach by utilizing the scheme for 
the Dicke model, for two applications. We have applied 
feedback to provide protection (against loss) of quantum 
states of interest for metrology (NOON and entangled 
coherent states). We have also demonstrated that a re¬ 
alistic feedback proposal including a time-delay can be 
used to alter the properties of the NESS potentially in¬ 
creasing the entanglement between the sub-systems and 
changing the stability of the steady state. 
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Appendix A: Markovian master equation with 
optical feedback 


In this appendix we provide the details of the deriva¬ 
tion of the master equation (21). For the feedback loop to 


be effective, the driven cavity should respond much faster 
than the source cavity. We thus work in a regime where 
7 d ^ 1 and the state of the driven cavity follows adiabat- 
ically the evolution of the source cavity. The large decay 
rate and the coupling to a zero temperature reservoir 
implies that driven cavity may only be weakly excited, 
which justifies the ansatz state ( [2Q| ). With such a joint 
state of two coupled cavities one obtains the following 
equations of motion for the source cavity 


7 ^ 

+ + IdPll 


poo — ^nAD,P00 

+\/77d(^Pio + Pio^^) + + PlO^) 


pio — 


7 - 7d - 

^nAD,PlO + -^f^aPlO - ^^10 


+\/TM(a/5ii - apoo + 

(v^Q^/ 020 - a/3oo + Piia) 


Pii = 


H, 


nAD, Pll 


+ -^^aPll ~ (^^d + 7d)Pll 


-Vrm(a/3io + PW^O - P^^^^iaPw + Pw^^) 


P 20 = — 


^nAD, P 20 


+ -j^^ap20 — 27dP20 


-y277dapio - p ''^^^ apw. 


dynamics 

P = Poo + Pii 


-^nAD 5 P 



(A4) 


where 7 eff = 7 (l + /i (2 + /i)) is the effective damping rate 
of the field mode in the source cavity. 


Appendix B: Uhlmann fidelity in the RWA regime 

A conspicuous feature of Fig. [^(especially in (c)) is the 
almost constant evolution of the fidelity for the states 
with N = 1, 2, and the oscillatory structure for the 
A = 3, 4 states. This implies that are station¬ 
ary states while seems to be formed from two sta¬ 

tionary states. In this appendix we will explain how this 
comes about given the initial states and the Hamiltonian 
( [ 1 ^ . As pointed out in the main text, for the figure a 
small coupling has been used (20 times smaller than the 
bare frequencies) which means that imposing the RWA 
is justified, i.e. we let A 2 = 0 from now on. We have 
numerically verified that the RWA is applicable for the 
corresponding parameters. 

Now, when A 2 = 0 the Hamiltonian can be readily di¬ 
agonalized by defining two new bosonic operators; ‘even’ 
= (d + h)l\/2 and ‘odd’ Z- = {a — h)/\/2. These 
obey the regular boson commutation relations and mu¬ 
tually commute. In particular the ‘odd’ Fock states 
|n_) = | 0 )/y^nin have shown to be important for 

adiabatic passage in multimode cavities m- Expressed 
with the new rotated operators, the Hamiltonian is diag¬ 
onal 


The state of the source cavity is of interest to us and 
it can be extracted as p = Tr^W = poo + Pii- From the 
above equations, the off-diagonal elements pio and p 2 o 
can be adiabatically eliminated by slaving them to the 
diagonal elements poo and pn. Setting p 2 o =0 we obtain 
to leading order in I / y/(7d / 7 ) 

P 20 =-y . , . (7 + l)apiQ. (AI) 

V2(7d/7) ^ 

Inserting the above expression in the steady state equa¬ 
tion of pio we obtain 

==(—Pun — —CLpoo + npii — apoo)- (A2) 

h) ^ ^ 

Using the above steady state solution for the off-diagonal 
element yields the following master equation for the den¬ 
sity matrix of the source cavity alone 

p = Poo + Pll = —p] + — >CaP — iUdPii 

+P7[>Ca(poO ~ Pll) + ~>^aP 00 + Pll] •(A3) 

When 7 d ^ I, pn ^ 0(0) and p ~ poo we arrive at 
the following master equation approximating the source 



Hz — (c(j -f Ai)z^z+ + (cj — \i)z^_Z—. (BI) 

It follows that the eigenstates are |n+)|n_) = 

y^n^JnZT, with corresponding eigenen- 
ergies 5n+n_ = ^ (^+ + n_) -h Ai (n+ - n_). 

For an N particle NOON state we have 


|^7 = ^(at^- 6 tivyo)JO),. (B2) 

Thus, we notice that |TL) = ii|0)a|0)6 is indeed a sta¬ 
tionary state with its time evolution given by 




(B3) 


The time evolved state (B3) is of course strictly only cor¬ 


rect within the RWA, and the exact time evolved state 
should be given by evolution under the full Hamiltonian; 
l^exact(^)) = Numerically we find the er¬ 

ror arriving from neglecting the counter rotating terms 
= Inexact ~ ^rwaI < 0.003 for all t. Similarly, we find 
that the two particle NOON state |T?_) = z\_z^\0)a\0)b 
is also a stationary state with the time evolution 


(B4) 
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In this case the numerically estimated error S 2 = | ^exact ~ 
^approxI < 0.004 for all t. Interestingly, the alternative 
NOON state \^\) = (|2)a|0)6 + \0)\a2)b)/^/2 is not a 
stationary state. 

Turning now to the N = 3 case. It is easy to show 
that there is no A/" = 3 NOON state (most general form 
13)a 1 0)5 + e^^\0)a\3)b) that can be a stationary state. For 
the present state the time evolved state takes the form 

p—iSojt 

2v3! 

(B5) 

It is a straightforward exercise to show that |T^(t)) 
evolves as a time dependent mixture of the orthogonal 
states ) and l^'^th) = (|2)a|l)b - |l)a|2)b)/V2. But, 
never in the course of time evolution the coefficient of 
|T^) vanishes completely. This is the reason why the 


fidelity 4> for = 3 in Fig. (c) never drops close to 
zero. 

The situation is however somewhat different for 
|Tl(t)) whose time evolution is 

v4! 

(B6) 

Even though, like in the N = 3 situation, one cannot 
find a general NOON state being a stationary state (in 
fact it is only possible for N = 1, 2), the above time 
evolved state evolves into a state completely orthogonal 
to |Ti). For instance, when 2Ait = (2m + I)7r/2 one ob- 

tains |«'i(t)) - = (|3)a|l)b - \l)a\^)b)/V2. This 

is where the fidelity 4> for |Tl) drops close to zero as 
shown in Fig. (c). 
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